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The Hubble Deep Fields continue to be a valuable resource for studying the 
distant universe, particularly at z > 2 where their comoving volume becomes 
large enough to encompass several hundred L* galaxies or their progenitors. Here I 
present recent results from a near-infrared imaging survey of the HDF-North with 
NICMOS, which provides structural and photometric information in the optical 
rest frame (AAo4000-5500A) for hundreds of ordinary galaxies at 2 < 2: < 3, and 
which offers the means to search for still more distant objects at z ^ 5. Lyman 
break galaxies at 2 < z < 3 are compact and often irregular in the NICMOS images; 
ordinary Hubble sequence spirals and ellipticals seem to be largely absent at these 
redshifts, and apparently reached maturity at 1 < z < 2. The Lyman break galaxies 
have UV-optical spectral energy distributions like those of local starburst galaxies. 
Population synthesis models suggest typical ages ^fewxlO® years and moderate 
UV extinction (~ 1.2 mag at 1700A), but the constraints are fairly weak and there 
may be considerable variety. Considering a near-IR selected galaxy sample, there 
is little evidence for a significant number of galaxies at z '--^ 3 that have been missed 
by UV-based Lyman break selection. Using the well-characterized z ~ 3 galaxy 
population as a point of reference, I consider Lyman break galaxy candidates at 
4.5 < z < 9, as well as one remarkable object which might (or might not) be at 
z > 12. The space density of UV-bright galaxies in the HDP appears to thin out 
toward larger redshifts, although surface brightness selection effects may play an 
important role. 
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1. Introduction 

(a) The first galaxies? 

The past five years have seen remarkable breakthroughs in our ability to identify 
and systematically study ordinary galaxies at very large redshifts, not just as iso- 
lated case studies, but en masse as a galaxy population. To date, nearly 1000 galaxies 
have been spectroscopically confirmed at z > 2, mostly identified via broad-band 
color selection techniques (Steidel et al. 1996, i.e., the Lyman break galaxies, or 
LBGs), but with other important and complementary methods also coming into 
play (sub-mm and radio surveys, emission line searches, QSO absorption systems, 
etc.). A broad census of galaxy properties at z ~ 3 now seems within reach, cov- 
ering star formation rates (SFRs), dust content, morphologies, spatial clustering, 
and perhaps even chemical abundances and internal kinematics. 
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Although I have retained my assigned title, 'The First Galaxies...,' it is far from 
clear that we know when, where, or how to find the 'first' galaxies. The 2; ~ 3 LBGs 
may or may not be the first major wave of galaxy formation. If anything, current 
data favors a roughly constant global SFR (as traced by cosmic UV luminosity 
density, at least) from 2 < z < A, with no certain evidence for a decline at higher 
rcdshifts (Stcidcl et al. 1999). The sub mm population detected by SCUBA (cf. 
Cowie, this volume) may or may not represent the bulk of early star formation, and 
the upper redshift bound to SCUBA sources remains unknown. The reionization of 
the intergalactic medium at z > 5 and the presence of metals in the Ly a. forest 
at 2; ~ 3 point toward earlier epochs of star and galaxy formation, at least in trace 
amounts. A handful of galaxies now have plausible spectroscopic confirmation at 
> 5, but too few for any systematic census. 

For this reason, I cannot promise to live up to my title: I do not know what 
the first galaxies are or what they look like. Given this ignorance, I will focus 
the first part of this article on the structure and stellar populations of the most 
distant well studied galaxies, the 2 ~ 3 Lyman break objects. This is not meant as a 
comprehensive review, but will instead concentrate on new imaging and photometric 
data from NICMOS on HST that extend our knowledge of LEG properties to 
the optical rest frame. In his contribution to this volume. Max Pettini provides a 
complementary discussion of recent efforts to measure chemical abundances and 
internal kinematics for these same galaxies. In the second part of my article, I 
will describe efforts to extend Lyman break color selection to still larger redshifts, 
approaching (or perhaps even exceeding) z ~ 10. In this way I hope to at least 
provide a look into the epoch when 'the first galaxies' might plausibly have been 
formed, and to catalog what we can find right now, in the pre-iVG^Tera, given the 
best available survey data. 

(6) Infrared observations of the Hubble Deep Field 

For the past five years, the Hubble Deep Fields (HDFs) have provided the most 
exquisitely deep, high angular resolution optical census of the distant universe. It 
is important (if somewhat pedantic) to consider what an HDF is actually good for. 
One WFPC2 field covers 5 arcmin^, and probes a very small co-moving volume at 
z <1, enough to hold only ~ 12-30 L* galaxies, depending on the cosmology. Given 
small number statistics and concerns about clustering, the central HDF is therefore 
not the best place study massive galaxies in the 'low' redshift universe, despite the 
fact that most of cosmic time and most bright galaxies with spectroscopic redshifts 
are at 2; < 1. There is far more volume at high redshift: 10.5 to 40 x more at 
2 < z < 10 than at z < 1 for plausible cosmologies, room enough for several 
hundred L* galaxies or their progenitors. 

At 2 > 1, the optical light emitted from galaxies shifts into the near-infrared. 
Thus in order to compare z > 2 galaxies to their local counterparts, and to search 
for still more distant objects at z ^ 5, it is important to extend the wavelength 
baseline. The HDF-North was observed in the near-IR from the ground in several 
different programs (e.g.. Hogg et al. 1997; Barger et al. 1998; and our own KPNO 
4m JHKg imaging, cf. Dickinson 1998). The depth and angular resolution (typically 
1") of these data are a poor match to that of the optical WFPC2 HDF images. 
Two programs therefore targeted the HDF-North with NICMOS on board HST, 
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providing much deeper images with high angular resolution. The NICMOS GTOs 
(Thompson et al. 1999) imaged one NICMOS Camera 3 field (~ 51" x 51") for 49 
orbits each at FllOW (1.1/im) and F160W (1.6/im). Our own program mosaiccd 
the complete HDF with a mean exposure time of 12600s per filter in FllOW and 
F160W. Sensitivity varies over the field of view, but the average depth is AB « 26.1 
at S/N = 10 in an 0"7 diameter aperture. The drizzled PSF has FWHM = 0"22, 
primarily limited by the N1C3 pixel scale. Because most galaxies have spectral 
energy distributions (SEDs) which brighten (in f^, units) at redder wavelengths, 
our images detect roughly half of the galaxies from the WFPC2 HDF, despite their 
short exposure times. We have also reanalyzed our KPNO Kg images to optimally 
extract photometry matched to the WFPC2+NICM0S data. These data are not as 
deep as one would like, which is unfortunate because they provide the only access 
to rest frame optical wavelengths for objects at 3 < z < 4.4, but they are the best 
presently available. 

Thanks to the dedicated efforts of the observers, a remarkably high density of 
spectroscopic redshifts is available in the HDF-N: ~ 150 galaxies (plus a few stars) 
in the central WFPC2+NICM0S field alone, with 33 objects at 2 < z < 5.6. Tak- 
ing advantage of the high quality photometric data, many investigators have used 
multicolor selection (e.g. the two color Lyman break criteria of Steidel et al. 1996, 
Madau et al. 1996, and others) or photometric redshifts (e.g., Fernandez-Soto et 
al. 1999) to identify high rcdshift galaxy candidates. There arc advantages and 
drawbacks to both approaches, but both have demonstrated remarkable successes. 
In this discussion, I will make use of both methods. For the photometric redshifts, 
I will use fits to our 7 band HDF photometry by Budavari et al. (2000), whose 
'adaptive template' method is a modification of an otherwise straightforward spec- 
tral template fitting scheme. I will use AB magnitudes here throughout, and notate 
the six WFPC2 + NICMOS bandpasses by U:mh, B450, Veoe, hu, Jiw and Hieo- 
Unless stated otherwise, I will assume a cosmology with = 0.3, = 0.7, and 
Ho = 70 kms-^Mpc-^ 

2. Galaxy morphologies at 2 < z < 3 

The NICMOS Hiqq images sample rest frame wavelengths in the optical V to B 
bands from z = 2 to 2.8. This upper bound is about the midpoint of the redshift 
range where C/300-dropout Lyman break selection in the HDF is most efficient. 
Therefore, at these redshifts we may use the NICMOS data to study the mor- 
phologies of LBGs at wavelengths where long-lived stars, if they are present, may 
dominate the light from the galaxy, and where dust obscuration should play a sig- 
nificantly lesser role than it does in the ultraviolet. For LBGs at z > 3, the NICMOS 
iJieo bandpass slips into the rest frame ultraviolet, and the NICMOS images once 
again tell us more about the distribution of star formation within galaxies than 
about that of their stellar mass. 

Figure 1 compares rest frame UV and optical images of a set of HDF LBGs 
with spectroscopic redshifts 2 < 2: < 3. The NICMOS data have somewhat poorer 
angular resolution (0"22 compared to 0"14 for WFPC2), but otherwise the most 
striking thing is the broad similarity of the UV and optical morphologies. Giavalisco 
et al. (1996) and Lowenthal et al. (1997) have emphasized the very small sizes of 
LBGs in WFPC2 images, and the same is true at NICMOS wavelengths. Accounting 
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Figure 1. Ultraviolet and optical rest frame images of HDF Lyman break galeixies with 
spectroscopic rcdshifts 2 < z < 3. For each object, the left hand image is interpolated 
between WFPC2 passbands to Aol700A, while the right panel interpolates the NICMOS 
images to the rest frame S-band (4300A). Galaxies at 2.8 < < 3 use the NICMOS i?i6o 
images without extrapolation. Boxes are 4" x 4", or 32/i,y(/ kpc on a side at z = 2.5. Here, 
the WFPC2 images have not been convolved to match the NICMOS PSF. 



for PSF differences, the half-light radii of the galaxies measured from the NICMOS 
images are the same as (or in some cases, slightly smaller than) those from the 
WFPC2 data. Notable morphological differences are seen in only a few cases, e.g., 
2-585.1 {z = 2.008) and 4-52 {z = 2.931), two of the largest LBGs in the HDF. 
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Each has a region of diffuse emission in the WFPC2 images which lights up in the 
NICMOS i?i6o data, with quite red colors. It is not clear whether this is due to the 
presence of dust, older and redder starlight, or possibly strong line emission.|t| 

In general there is no evidence that the UV-bright regions seen by WFPC2 are 
just star-forming 'fragments' embedded within some larger, mature host galaxy. 
On the whole, one or more of the following appears to be true: (1) the stars which 
dominate the light at AAol200-1800A also dominate at AAo4000-5500A; (2) if com- 
ponents with substantially different ages and colors are present, then they are fairly 
well mixed, spatially; (3) dust extinction does not play a dominant role in shaping 
the morphologies of LBGs at these wavelengths. All of this should be taken modulo 
the important caveat that at FWHM = 0"22 resolution, most of the LBGs are 
not exquisitely resolved: they are typically very small, having only a few resolution 
elements within their isophotally detectable areas, and thus many details are surely 
lost except perhaps for the few, largest objects. 

It is striking that among the NICMOS images of 100 HDF LBGs (those with 
and without spectroscopic redshifts), virtually none resembles a 'classical' Hubble 
sequence spiral galaxy. Giant disk galaxies are found in the HDF and elsewhere out 
to at least z « 1.3 (see, e.g., examples from NICMOS in Bunker 1999 and Dickinson 
2000), and various studies have found that their structural properties and comoving 
abundances have not changed dramatically since z « 1 (Lilly et al. 1998; Simard 
et al. 1999). At z > 2, however, there appear to be no objects with immediately 
recognizable spiral structure, no bulges surrounded by symmetric, diffuse disks, nor 
even good candidates for thin, edge-on disks (some LBGs are fairly elongated, but 
none really could be mistaken for an edge-on spiral). Among the objects shown in 
figure 1, perhaps 2-585.1 (z = 2.008) appears closest to showing spiral structure 
at rest frame optical wavelengths, but this requires some imagination. Given the 
evidence from NICMOS, the absence of classical spiral morphologies cannot be 
attributed solely to rest frame wavelength effects. We have carried out simulations 
using optical rest frame images of Virgo cluster spirals and of HDF disk galaxies 
at z < 1, artificially inserted into the NICMOS images at high rcdshift with the 
appropriate surface brightness dimming and PSF convolution. Although low surface 
brightnesses and limited angular resolution wipe out many of the details of spiral 
structure, giant [L > L*) disk galaxies should be detectable and recognizable at 
z > 2 even with no luminosity or surface brightness evolution. This is not to say 
that LBGs cannot be disk galaxies: as already noted, the small angular sizes of 
most LBGs preclude detailed resolution, and some of these objects could well be 
small disks. Indeed, Giavalisco et al. (1996) noted exponential surface brightness 
profiles among some LBGs. 

Giavalisco et al. also found that some LBGs have R^^^ law profiles, although 
it is not clear whether to interpret this as indicating a relation between LBGs and 
present-day ellipticals and bulges. It is worth noting, however, that there are few 
candidates for intrinsically red elliptical galaxies at z > 2 in the HDF, i.e., objects 
which would be recognized as elliptical galaxies today both by morphology and by 
the characteristic color signature of an older stellar population. In fact, considering 
objects with spectroscopic or photometric redshifts in the range 2 < z < 4, where 

t 4-52 is one of the few z > 2 galaxies detected at 8.5 GHz by Richards et al. 1998. The radio 
centroid is coincident with the diffuse, red, IR-bright region in the galaxy. 
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the NICMOS+KPNO infrared data still reach the optical rest frame, there are only 
a handful of galaxies with rest frame colors redder than that of a present-day Scd 
spiral (i.e., an actively star forming galaxy), even in an infrared-selected sample 
which should have no bias against such objects. I return to this point in §4 below. 
Nearly every HDF galaxy at z > 2 which is detectable at 1.6/im appears to be 
forming stars, and most quite vigorously. We do find apparently red, dead ellipticals 
in the HDF out to (photometric) redshifts z 1.8 (cf. Dickinson 2000; Stanford 
et al. in prep.), but only one (marginally) viable candidate for a red elliptical at 
higher redshift. This is the so-called ' J-dropout' object HDFN-JDl (see §6), whose 
colors might be matched by those of a maximally old elliptical galaxy at 3 < z < 4 
(Dickinson et al. 2000) .| 

Overall, it seems that the maturation of the giant spiral and elliptical galaxies 
took place at z < 2. Even with extremely deep, high angular resolution infrared 
images like the HDF/NICMOS data, we find few (if any) mature spirals or ellipticals 
(or candidates for such) at higher redshift. By z — 1, many (if not necessarily all) 
large spirals and red giant ellipticals were already in place, pointing to the redshift 
range 1 < z < 2 as an important 'golden age' for the formation of the Hubble 
sequence. 

3. Galaxy colors at 2 < z < 3.5 

Figure 2 shows a photometric compendium of spectroscopically confirmed HDF 
galaxies at 2 < z < 3.5, all shifted to the rest frame and normalized to a common 
flux density at Ao2000A. The shaded envelope is defined by local UV-to-optical 
starburst galaxy spectral templates from Kinney et al. (1996), which span a broad 
range in optical/UV extinction. The HDF LBGs fall comfortably within the range of 
SED shapes defined by the local starbursts. They have relatively blue (but usually 
not flat spectrum) UV continua, with a flux increase and spectral inflection around 
the Balmer/4000A break region that indicates the presence of older (A and later) 
stars which apparently contribute a significant fraction of the rest frame optical 
light. The fact that the UV continuum slope for LBGs is nearly always redder than 
flat spectrum has generally been interpreted as an indication of dust extinction 
(e.g., Meurer et al. 1997, 1999; Dickinson 1998; Pettini et al. 1998), although for 
some objects it might also result from an aging stellar population with declining or 
inactive star formation. Considered individually, the large majority of HDF LBGs 
are reasonably well fit by the Kinney et al. starbursts with modest reddening (0 < 
E{B — V) < 0.21). Few approach the more heavily reddened templates. This may 
be true by definition/selection, of course, since we are considering the brighter 
objects for which redshifts were successfully measured, and which were selected for 
spectroscopy by their UV colors (but see §4 below). 

With photometry spanning the UV-to-optical rest frame it becomes interest- 
ing to compare the LBG photometry to population synthesis models to look for 
constraints on galaxy ages, reddening, and star formation histories. Very roughly, 
if the UV spectral slope provides a measure of extinction (modulo an assumed red- 
dening law), then the UV-optical flux ratio, and particularly the amplitude of any 
inflection around the Balmer/4000A break region, may help constrain the past star 

f In fact, HDFN-JDl is probably too red for an old stellar population at high redshift without 
invoking dust, an unusual IMF, or an unfashionable cosmological model. 
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Figure 2. Photometry for 27 HDF galaxies with spectroscopic redshifts 2 < z < 3.5, shifted 
to the rest frame and normalized at 2000 A. The shaded region spans the range of empirical 
starburst templates from Kinney et al. (1996). The starburst SED sequence is primarily 
defined by reddening up to E{B — V) < 0.7; the blue envelope is set by NGC1705. 

formation history, particularly the ratio of older stars to ongoing star formation. 
This was first done by Sawicki & Yee (1998) using ground-based JHKg photometry 
for LBGs in the HDF. Their results favored young ages (median value ^ 25 Myr) 
and fairly heavy reddening (typical E{B — V) fa 0.28, or 3 mag extinction at 1600A 
assuming Calzetti 1997 starburst dust attenuation). We have carried out a similar 
exercise using the NICMOS data and i^-band fluxes rederived from the KPNO 
data using a technique (much like that of Fernandez-Soto et al. 1999) which prop- 
erly matches photometry from images with very different angular resolutions. A 
complete presentation will be given in Papovich et al. (in prep.) and is beyond the 
scope of the present discussion, but I summarize some important points here. 

Even with precise NICMOS photometry, constraints on ages and reddening are 
quite loose. This is in part because of the usual degeneracies in fitting models to 
broad band colors (age vs. metallicity vs. extinction), but also because the available 
photometry simply does not reach long enough rest frame wavelengths. At z > 3, 
only the \owei-S/N ground-based data extend redward of the Balmer/4000A 
break region; for this reason, here I consider only galaxies at 2 < z < 3. For 
each galaxy, we may define confidence intervals in the multidimensional space of 
the various fitting parameters such as age, SFR e-folding timescale, reddening, 
metallicity, and IMF. The color degeneracies allow fairly broad ranges in acceptable 
parameters for each object, and the variety among the galaxy SEDs (figure 2) 
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Figure 3. Combined confidence intervals (68%, 95%, 99.7%) for population syntfiesis 
model fits to a sample of 16 HDF Lyman break galaxies with spectroscopic redshifts 
2 < z < 2.95, from Papovich et al. (in prep.). A wide range of Bruzual & Chariot (1996) 
models has been considered with varying age (t), star formation histories (exponential 
timescale r), dust attenuation {E{B — V), here assuming the Calzetti 1997 starburst law), 
metallicity, and IMF. Here, results for the 16 objects have been averaged, weighted by 
their individual probability distributions, to show the ensemble likelihood distributions 
for E{B — V), t and r. Best-fitting values for each object are marked by points in the 
left panel. Dependences on metallicity and IMF have been collapsed, but are fairly weak. 
The cross-hairs at left indicate the most-favored values for E{B — V) vs. t. The dashed 
line at right marks t = t. Models observed after the bulk of their star formation has been 
completed (i.e. t S> t), are disfavored except at young ages (< few x 10* years) and short 
timescales, when the UV continuum can persist throughout the main sequence lifetimes 
of B stars. 

scatters the best-fitting parameters throughout a range of values. Nevertheless, 
taken as an ensemble, certain regions of model parameter space are preferred. 

Figure 3 shows a composite distribution for 16 galaxies with z < 3, where 
model fits for each object have been averaged, weighted by their likelihoods in 
the multi-parameter space. The contours thus indicate a distribution of likely pa- 
rameter values for the ensemble of UV-selected LBGs. The best-fitting values of 
E{B — V) vs. age for individual objects are marked by dots. Some galaxies can 
be fit reasonably well by parameter values falling toward the outer contours of the 
ensemble distribution, but the majority occupy the higher confidence regions. The 
most favored age range spans 0.03-1 Gyr, with extinction < E{B — V) < 0.25. 
The extinction values agree well with the comparison to the Kinney et al. starburst 
templates (figure 2), which is not unexpected given that the Calzetti attenuation 
law is derived in part from the same UV spectral data on local starburst galaxies. 
The range of likely age and extinction values becomes slightly smaller and can shift 
somewhat if restrictions on metallicity or IMF are adopted, although in general the 
broad band photometry offers little constraint on these parameters. Younger ages 
and larger extinction are allowed and even favored for some objects, but in general 
the most likely values are somewhat older and less reddened than those found by 
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Sawicki & Yee. This may be due to more precise photometry from NICMOS, or to 
better control of the relative optical-IR colors (particularly in the ii'-band) from 
our photometric method. 

It is not clear whether the apparent anticorrelation between age and extinction 
in figure 3 is significant. We expect degeneracy between the age and extinction 
values fit to individual objects, but the overall trend seen in the population may 
be greater than would be expected from the fitting uncertainties alone. The 'most 
favored' extinction value, E{B — V)w 0.12, corresponds to A(1700A) f« 1.2 mag, 
or a factor of ~ 3. The net UV extinction for the sample (and hence the correction 
to any derived global star formation rate) would be larger, however, driven by the 
objects with the greatest reddening. 

A characteristic timcscalc for Lyman break galaxies can be defined from their 
sizes (median half-light radii fa 2.2 kpc for the HDF LBGs) and typical velocity 
dispersions (~ 80 km s~^; see Pettini, this volume). This yields tc ~ 25 Myr. If the 
UV light is due to ongoing star formation, we would not expect SFR lifetimes <C t^, 
and indeed this is roughly the lower bound of the best-fit model age range, with 
the most favored value being ~ 135 Myr, or ~ 5tc- This range is not dissimilar to 
that estimated for star formation in galactic-scale starburst events (e.g., ULIRGs) 
locally. 



4. Rest frame UV selection: what do we miss? 

Using an infrared-selected catalog, wc may ask what galaxies might be missed 
altogether by Lyman break color selection keyed to the rest frame UV light. In 
particular, one might expect some red high redshift galaxies, either because they 
arc not actively forming stars or because of extinction, that would 'drop out' of 
the dropout samples. Here I restrict my analysis to Hiqq < 26, where we believe 
our catalogs are highly complete, uncontaminated by spurious sources, and where 
the NICMOS photometry has S/N > 10. At z = 2.75, iJujo < 26 corresponds 
to rest frame Mb < —19.36 for the adopted cosmology, or ~ 1 mag fainter than 
present-day Lg. The typical LBG at Hieo ~ 26 has Veoe ~ 27, the practical limit 
for HDF UsoQ dropout selection using standard 2-color criteria, but red galaxies 
with similar rest frame optical luminosities might be fainter or absent in the UV. 

At ^160 < 26 there is only one object which is undetected with S/N < 2 in 
Veoe or /814 (both, in this case): this is the 'J-dropout' HDFN- JDl (see also §§2 
and 6). In fact, this is the only NICMOS- selected object with -ffieo < 26 and 
S/N{Isi4) < 6.5. Two other objects have S/N{Vgo6) < 3; both are 2; > 5 Veoe- 
dropout candidates identified by Lanzetta et al. (1996) and Fernandez Soto et al. 
(1999), one of which (3-951) was spectroscopically confirmed at z = 5.33 (Spinrad et 
al. 1999). Thus the only possible candidate for a NICMOS-selected, 'UV-invisible' 
galaxy at ;j - 3 is HDFN-JDl. 

Next let us consider UV-bright objects which might nevertheless have been 
missed by the LBG color criteria, using the 7-band photometric redshift estimates 
for all galaxies. In principle, these may identify plausible candidates at 2 < z < 3.5 
that otherwise fall outside a given set of UV color criteria, as long as their intrinsic 
SEDs are 'recognizably similar' to those of galaxies at lower redshift which define 
the templates used for the phot-2; fitting. 
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Figure 4. Photometric and spectroscopic redshift distributions for infrared-selected HDF 

galaxies {Hieo < 26) at 1.5 < z < 4. The open histogram shows the Zphot distribution 
from Budavari et al. (2000). The hatched histogram indicates galaxies which obey the 
Lyman break criteria defined in the text, while the filled histogram shows the available 
spectroscopic redshifts. 



Figure 4 shows the Zphoi distributions for all HDF galaxies with 21 < iJieo < 26, 
and for those that meet the C/300 dropout criteria {Uzoo — B450) > {B45Q — Veoe) + 1 
and B450 — Veoe < 1-2 from Dickinson (1998). There are 43 objects with 2 < 
Zphot < 3.5 which do not meet the LBG color criteria. However, nearly all are 
at 2 < Zphot < 2.2 or 3.2 < Zphot < 3.5, and lie just outside the color selection 
boundaries defined here: either slightly too blue at low z, or slightly too red at high 
z. This is expected: the selection efficiency of the 2 -color method is not uniform 
with redshift, and falls off at the extremes of the range for which it is optimized 
(cf. Steidel et al. 1999). Only 7 'missed' objects fall at intermediate photometric 
redshifts, 2.5 < Zphot < 3.1, and most of these are also just outside the color 
selection box. Some are quite interesting, including a /zJy radio source with very 
red Jiio — Hieo colors which may be a dusty starburst or a fading post-starburst 
galaxy at z ^ 2.6. Others may be scattered out of the box by photometry errors 
(especially in t/300), or might not be at the indicated Zphot- But all are well detected 
in the optical HDF. 

Overall, there is no evidence for a substantial population (by number) of galax- 
ies at 2 < z < 3.5 that are missed by UV Lyman break color selection but which 
are detectable in the near-infrared. If there are energetically important but highly 
obscured galaxies at these redshifts like those detected by SCUBA, then they are 
either also detectable with optical imaging data, or they are so heavily enshrouded 
that even NICMOS cannot easily see them. For the five sub-mm sources detected 
in the HDF by Hughes et al. (1998), our NICMOS images do not reveal any new 
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Figure 5. Apparent magnitude for an L* (AqITOOA at 2: « 3; Steidel et al. 1999) Lyman 
break galaxy redshifted without evolution. No fc-correction is included: instead, the de- 
tection bandpass is taken to be fixed to Aol700A and changes with redshift (see top axis 
labels). An L* LBG should be detectable in the NICMOS data (i/ieo < 26.5) out to 
2 ~ 10 for spatially flat cosmologies with A < 0.8, and to z k. 7.5 for an open universe. 

counterparts previously undetected by WFPC2, nor do any of the candidate iden- 
tifications have particularly unusual optical-IR colors. 

5. Galaxies at 4.5 ^ z ^ Q 

The successes of color selection techniques at 2 < z < 4.5 make it tempting to 
extend the methods to higher redshifts, i.e., to search for V- or /-dropouts. Doing 
so properly requires deep near-infrared data to provide at least one color longward 
of the redshifted 912A and 1216A breaks. Indeed, galaxies at z > 6.5 should have 
virtually no detectable optical flux. Lanzetta et al. (1996) and Fernandez-Soto et 
al. (1999) identified candidate z > 5 HDF galaxies from Veoe — hu colors supple- 
mented by infrared limits from the KPNO data; two of these have subsequently 
been confirmed via spectroscopy (Weymann et al. 1999; Spinrad et al. 1999). With 
NICMOS we can extend this to fainter limits and larger redshifts: an L* Lyman 
break galaxy (i.e., L* in the rest frame UV at z = 3) should be detectable in the 
NICMOS images with -ffieo < 26.5 out to z = 10 for spatially flat cosmologies with 
A < 0.8, and out to z « 7.5 for an ^Im — 0.2 open universe (see figure 5). 

Here we test the null hypothesis that the galaxy population at z 3> 3 is similar 
to that of the J/soo-dropout LBGs at z ~ 3, whose observed characteristics are by 
now reasonably well known even if their intrinsic properties, such as dust content, 
star formation rate, mass, etc. are the subject of continued debate. In particular, we 
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Figure 6. Two-color diagrams for combinations of adjacent HDF WFPC2+NICM0S fil- 
ters. Symbol size scales with magnitude; known and suspected stars and a few bright 
galaxies have been excluded for clarity. Triangles mark Icr lower color limits; filled sym- 
bols indicate galaxies with spectroscopic redshifts. The lines show the nominal color-vs-z 
tracks for various unevolving galaxy SEDs, and the shaded region indicates the approx- 
imate color range expected for galaxies in redshift ranges appropriate to each color pair 
{UBV: 2-3.5; BVI: 3.5-4.5; VIJ: 4.5-6; IJH: 6-8.5). The color selection boxes used for 
the comparison with Monte Carlo simulations are indicated. 



adopt the rest frame UV luminosity function and UV spectral slope (i.e., intrinsic 
color) distribution for LBGs at (z) « 3 derived in Steidel et al. (1999), and use 
this to predict what should be seen in color-color diagrams ij the same population 
were present at higher redshifts. We do this via Monte Carlo simulations, including 
realistic errors for the HDF WFPC2-I-NICM0S photometry, comparing the number 
of high-z objects that are predicted to fall in some specified color-color box to the 
actual number of similar objects found in the HDF catalogs. 

Figure 6 shows a series of 2-color diagrams for the HDF, each using combinations 
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JjjQ magnitude H,gg magnitude 

Figure 7. Cumulative number counts of objects satisfying the high redshift color criteria 
shown in figure 6. The irregular histograms are from the HDF data, and the smoother 
histograms are predictions from the Monte Carlo simulations which assume that the galaxy 
population at high redshift is like that observed at 2: ~ 3. Models were computed for 
two A = cosmologies: fiji/ — 1.0 and 0.2 (the open cosmology is the lower model in 
the VIJ and IJH plots). The vertical dashed line indicates a rough confidence limit in 
magnitude for each plot, below which the data should become significantly incomplete 
and/or contaminated by spurious sources. 

of three adjacent bandpasses, from UBV (i.e., z ~ 3 selection) through IJH (i.e., 
z ~ 7). In each case, I define somewhat arbitrary selection boxes based on the 
expected location of high redshift galaxies in color-color space (and also to avoid 
low redshift contaminants), then count the galaxies in those boxes and compare 
this to the 'no evolution' (NE) model predictions (figure 7). 

The [/300-dropout counts agree well with the models by construction, since 
the input luminosity function is partially based on HDF data. The i?45o-dropouts 
fall below the NE predictions. This is just the original Madau et al. (1996) result 
revisited: the HDF-N appears to have fewer galaxies at z ^ 4 than at z ~ 3. 
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Figure 8. UV luminosity at Aol700A vs. co-moving volume out to redshift z for galaxies 
with Hieo < 26.5. A constant horizontal density of objects implies a constant space density. 
Fluxes are measured through 12h^Q kpc diameter metric apertures at all redshifts. The 
right-hand axis indicates the corresponding 'raw' star formation rate derived from the UV 
luminosity without correction for dust, assuming a Salpeter IMF. The horizontal dotted 
line marks the luminosity of an L* LBG at z w 3 (Steidel et al. 1999). The dashed curve 
indicates m(Aol700) = 26.5 at the bandpass indicated by labels at top. Objects with 
spectroscopic redshifts arc indicated with filled symbols; photometric redshifts are used 
otherwise. Objects meeting 2-color criteria as U300, -B450, Veoe and /814 dropouts are coded 
as V) CHi and Oi respectively; all others are plotted with + and o. 



Steidel et al. (1999), who surveyed larger solid angles in several fields, suggest that 
the bright end of the 2; ~ 4 LF is actually compatible with that at 2; ~ 3. The HDF 
may just be an anomaly, indicating the importance of field to field fluctuations, or 
perhaps the faint end slope of the LF (to which the HDF number counts are quite 
sensitive) evolves with redshift. For the Veoe-dropouts, there are ~ 7 candidates 
with Jiio < 26.5 (including the two with spectroscopic confirmation), compared to 
a prediction of ~ 17. Careful inspection of their images and SEDs suggests that 
they are all very plausible 4.5 < 2; < 6 candidates. There are no 78i4-dropout 
candidates with iJieo < 26, and only two with -ffigo < 26.5, one of which is clearly 
detected at S450 and Veoe and thus is probably not at 2; > 6. The models predict 
9 to 13 objects to this magnitude limit. Some of the fainter objects may be real 
2 > 6 galaxies, but on visual inspection many are rather dubious, with very low 
S/N; only a few are persuasive to a skeptical eye. At i?ieo > 26.5 we are reaching 
or passing the useful depth limits of our NICMOS data for this purpose. 

This analysis is in qualitative agreement with one based on photometric redshift 
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estimates. Figure 8 plots rest frame 1700A luminosities of galaxies vs. redshift 
(spectroscopic when available, photometric otherwise) for an HDF sample limited 
to HiQQ < 26.5. The photometric redshifts are generally in good agreement with 
the simple 2-color selection illustrated in figure 6. There are a few objects with 
Zphot ~ 3-5 that 'fall in the gap' between the U^qq- and i?45o-dropout samples. At 
Zphot > 6 there is only partial overlap between the /814-dropout and photometric 
redshift samples, but this is understandable since most of the candidates are very 
faint with \ow-S/N photometry and poor photometric redshift constraints (i.e., 
relatively flat Zphot likelihood functions). The space density of the bright LBG 
candidates appears to thin at z > 4.5, and at z > 5.5 there are no candidates 
with UV luminosities greater than the characteristic L* at z = 3, despite abundant 
volume to house them if they were present with similar space densities.^ 

Ferguson (1998) and Lanzetta et al. (1999) have stressed the importance of 
cosmological surface brightness dimming when characterizing the galaxy population 
at > 3. This can affect the likelihood of detecting high redshift galaxies, as well 
as the fluxes measured with isophotally-based photometry schemes. Lanzetta et al. 
show that the global rate of star formation occurring in the regions with the highest 
UV surface brightness rises steeply with redshift, and argue that far more UV light 
may be present in 2; > 5 galaxies at fainter, unmeasured isophotal thresholds. We 
may partially address this by measuring fluxes and luminosities non-isophotally, 
e.g. by using apertures scaled by image moments or with fixed metric sizes. The 
photometry in figure 8 uses 12 kpc metric apertures, and thus is insensitive to 
surface brightness limits except as far as they affect galaxy detection in the first 
place: most faint object cataloging packages use isophotal detection thresholds. We 
have examined this by taking WFPC2 images of 2; ~ 2.5 LBGs, artificially shifting 
them to higher redshifts, and re-inserting them into the NICMOS data to assess 
their detectability. Galaxies like the L > L* LBGs at 2; ~ 3 should be detectable 
to at least z w 7 at the depth of our GO NICMOS images, and more easily in 
the HDF-South NICMOS field or the HDF-North GTO NICMOS image, each of 
which goes ^ 1 magnitude deeper than the data set discussed here. But a general 
census at such redshifts might indeed be woefully incomplete. 

Overall, the HDF data disfavor the null hypothesis that galaxies like the bright 
LBGs at 2 < z < 4.5 are present at z ^ 5 with similar space densities. The higher 
redshift galaxies are apparently either fainter, more rare, have lower surface bright- 
ness, or some combination thereof. At any rate, they are certainly more difficult to 
detect and study, at least in abundance, even with NICMOS. 

6. An object at z 12? 

Although the few /-dropout candidates in the HDF are very faint, paradoxically 
there is one comparatively bright ' J-dropout' object, shown in figure 9. Lanzetta et 
al. (1998) identified five possible sources in the KPNO Kg images of the HDF which 
were invisible in the WFPC2 data. Of these, four are undetected by NICMOS. One, 
however, which we call HDFN-JDl, has a robust 1.6/im detection (i/ieo ~ 25.2), 
but is at best only marginally {S/N < 2) detected at 1.1/im and shorter wave- 

f For an open universe, the higher Zpf^^t candidates would be more luminous and the upper 
envelope to their UV luminosity would be nearly flat, but the space densities would be even more 
sparse compared to those at 2 < 2: < 3.5. 
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U300 B450 V506 1814 J110 Hi 50 Kg 

Figure 9. HST and Keck images of HDFN-JDl at 0.3-2. 16^m. The field of view of each 
panel is 4" x 8". HDFN-JDl is identified by tick marks in the i/ieo and Ks panels. The 
Keck Ks image has been smoothed by a Gaussian with FWHM = 0"38. 



lengths. This optical 'detection,' if real, would be important, as it would probably 
exclude the most exotic hypothesis for this object, i.e., that it is a galaxy or QSO 
at z > 10. It will be difficult, however, to obtain much deeper optical data than the 
existing HDF WFPC2 images to provide a stricter limit. The red -ffieo — Kg color 
suggests z « 12.5 under the high redshift hypothesis, with the Lya forest partially 
suppressing the i?i6o flux. We obtained an iJ-band spectrogram of HDFN-JDl 
with CRSP at the KPNO 4m, and to our surprise detected a moderately convinc- 
ing emission line at 1.65/xm that could plausibly agree with Lya at z = 12.5. The 
line did not reproduce, however, in a subsequent reobservation at higher dispersion 
(see Dickinson et al. 2000 for the spectra and further discussion). 

If this is not a Lyman break object then it may be either heavily reddened 
and at arbitrary redshift (but most likely z > 2, given the colors), or possibly a 
maximally old elliptical galaxy at 3 < z < 4 (see §2). If it is really at z « 12.5, 
then it is either a galaxy whose unobscured star formation rate (computed from 
the UV luminosity) is several hundred Mq yr~^, or an AGN, perhaps one of the 
hypothesized population responsible for re-ionizing the universe. If so, however, 
such objects are rare at 2 < z < 13, with a space density several hundred times 
lower than that of present-day L* galaxies, and it is unlikely that most of today's 
galaxies began their life in such a way. 

Table 1 compares the UV luminosities of some confirmed and candidate galaxies 
at z > 4, compared to that of an L* Lyman break galaxy at z = 3. In some 
cases infrared data needed to measure fluxes at rest frame 1700A are not available, 
and estimates based on published photometric or spectroscopic fluxes redward of 
Lyman a have been used instead.^ 'Raw' star formation rates are computed from 
the UV luminosities assuming a Salpeter IMF and no dust obscuration. 

CDFa-Gl at z = 4.815 is the most luminous Lyman break galaxy yet identified 
(among both and G-dropouts) in our large, ground-based survey. The spec- 
troscopically verified 5 < z < 6 galaxies have luminosities that are mostly fairly 
typical of z w 3 LBGs, ranging from 0.66 to 1.5xi*. The z = 5.74 object from 
Hu, Cowie & McMahon (1999) has been described as 'extremely luminous,' but is 
actually quite typical for z w 3 to 4 LBGs, and somewhat fainter than HDF 3-951 
at z = 5.33. Chen et al. (1999) have identified a candidate galaxy at z — 6.68 from 

t L^{1700A) estimates for RDl are derived in part from a J-band measurement by Armus et 
al. (1998). For HDF 3-951 and 4-473 we use our own NICMOS Jno photometry. 
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Table 1. Properties of selected high redshift galaxies and galaxy candidates 

( Assumes Qm = 0.3, Ha = 0.7, h = 0.7, and that SFR = IMq year~^ produces 
Li,(1700A) = 9 X lO^'' erg s"^ Hz"^ Magnitudes are on the AB scale. ) 



name 


z 


m 


M 


L/L* 


'raw' 


SFR reference 






(UV at 


~ Aol700A) 


(Mq yr-i) 


L* LBG 


3.04 


24.5 


■21.1 


= 1 


12.8 


Stcidel et al. 1999 


CDFa-Gl 


4.82 


23.6 - 


■22.8 


4.7 


60 


Steidel et al. 1999 


RDl 


5.34 


25.5: - 


■21.0: 


0.9: 


12: 


Dey et al. 1998 


HDF 3-951 


5.33 


25.0 - 


■21.5 


1.5 


19.6 


Spinrad et al. 1998 


HDF 4-473 


5.60 


26.0 - 


■20.6 


0.66 


8.5 


Weymann et al. 1998 


HCMl 


5.74 


25.5: - 


■21.1: 


1.1: 


13.5 


Hu, Cowie & McMahon 1999 


A 


6.68 


23.9: - 


23: 


5.8: 


74: 


Chen et al. 1999 


HDFN-JDl 


12.5 ?? 


23.9 - 


■23.9 


13.7 


175 


Dickinson et al. 2000 



STIS slitlcss spectra. The available photometry is limited, but based on the spec- 
tral continuum flux density estimate, this object appears to be significantly more 
luminous than other known LBGs at 2 < 2; < 6, or than any < 2; < 10 candidates 
from the HDF/NICMOS sample. This is a remarkable result, if true, since the solid 
angle covered by the STIS field is < 1 arcmin^. On the other hand, if the J-dropout 
object HDFN-JDl were really aA, z = 12.5, it would be more luminous still, nearly 
3x brighter than CDFa-Gl at 2 = 4.82, with a 'raw' UV SFR = 175 Mq yr"!. 

7. Discussion 

Overall, the HDF/NICMOS data demonstrate both the promise and the challenges 
which lie ahead for finding and studying the 'first' galaxies. The rest-frame optical 
view of Lyman break galaxies presented in §§2 and 3 strongly suggests that the 
galaxy population at 2 < 2; < 3 had not yet achieved maturity. The giant, Hub- 
ble sequence spirals and ellipticals that dominate the high-mass end of the galaxy 
population today are not seen at z > 2. In a sample of HDF galaxies selected in 
the near-infrared, nearly all galaxies with spectroscopic or plausible photometric 
redshifts 2 < 2; < 3.5 are evidently forming stars quite rapidly and can also be 
identified via their emitted frame UV light. The evidence from SCUBA shows that 
there are occasional 'monsters' whose obscured star formation may be quite impor- 
tant to the global emissive energy budget from galaxies. The identification of these 
objects, relatively rare by number, remains an important dilemma. 

Broad band color selection has been the most successful means for identifying 
high redshift galaxies, but we seem to be pushing the limits of what can be accom- 
plished at 2 > 6 with present day capabilities. The NICMOS HDF images are the 
deepest near-IR data now available, and they do include plausible candidates for 
galaxies at 6 < 2 < 9, but they are relatively few, and most are quite probably 
too faint for spectroscopic confirmation. We probably should not expect to find 
galaxies much brighter than these candidates unless some of the 'first' galaxies were 
significantly more luminous than the boring old 'later' galaxies that we have now 
surveyed extensively at 2 « 3. This is not impossible of course: the Chen et al. 
object and HDFN-JDl are both possible (but unconfirmed) 2 > 6 candidates more 
luminous than any normal Lyman break galaxy at 2 < 5. Perhaps indeed there are 
very luminous, relatively unobscured proto-galaxies out there waiting to be found. 
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a hope that was once quite widespread, but which seems to have gradually faded 
in the modern era of 25th magnitude LBGs and optically invisible SCUBA sources. 
Perhaps it will still make a comeback... 

These few rather speculative candidates aside, the evidence from the HDF alone 
would suggest that the population of UV-bright LBGs may be thinning out at 
z > 5, at least for objects comparable to those at the bright end of the z « 3 
luminosity function. It should be remembered, however, that this was also the 
conclusion reached by Madau et al. (1996) at 2; > 3.5, a result that has since 
been challenged by larger surveys with extensive spectroscopy. It is undoubtedly 
dangerous to draw conclusions too strongly from one 5 arcmin^ field. However, 
extending this work to larger areas and more sightlines will be an expensive effort. 
Surface brightness dimming and limited solid angle coverage may limit our ability 
to see much more with NICMOS (assuming that it is successfully revived in 2001), 
and ground-based near-IR imaging may never go deep enough to detect any but 
the most luminous objects at z > 5. Wider fields imaged with the HST WFC3 
near-IR channel (coming circa 2004) may offer the best survey opportunity until 
NGST, but a substantial investment of observing time will be needed to survey 
adequate solid angles to sufficient depth. 

Alternatively, we may turn to other observing strategies, e.g., by taking advan- 
tage of gravitational lensing from foreground galaxy clusters to boost very distant 
objects to detectable magnitudes. Narrow band and blind multislit emission line 
searches are being carried out through airglow windows (e.g., at A9150A, corre- 
sponding to z » 6.5; cf. Crampton & Lilly 1999; Stockton 1999). Or perhaps con- 
certed efforts to identify SCUBA sources will indeed turn up objects at 2; ^ 5, 
where the advantage of the negative sub-mm fc-correction is enormous. 

These data are offering a first glimpse into the so-called 'dark ages,' and giving 
hope that there may be luminous things there to find and study. In some sense, we 
may not know that we've found the first galaxies until we can find no more beyond 
them. Holding to that standard will ensure that a more challenging (and hence, 
more rewarding) threshold of proof always lies ahead. 

I thank my HDF/NICMOS collaborators for their contributions and for permission 
to show results in advance of publication, particularly Ceisey Papovich, Tamas Budavari, 
Jean-Marc Dcltorn, Chris Hanley and Harry Ferguson. I also thank Chuck Stcidcl, Max 
Pettini, Mauro Giavalisco, Kurt Adelberger, and Alice Shapely for their long-standing 
and highly valued collaboration. I am grateful to the Royal Society for inviting me to 
this meeting and for supporting my travel and accommodations, and to the editors of the 
proceedings for their patience. This work was supported by NASA grant GO-07817.01- 
96A. 
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